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АЛГОРИТМ ОПИСАНИЯ ОВРАЖНЬ(ІХ ЦЕЛЕВЬІХ ФУНКЦИЙ 


Представлений алгоритм опису яружних цільових функцій, який пропонується використовувати для 
детальної візуалізації розповсюдження зон ураження сильнодіючих отруйних речовин у разі 
виникнення техногенної аварії. Така візуалізація є необхідною і можливою при застосуванні 
геоінформаційних технологій для оперативного реагування на подію. Розглянуті можливості 
апроксимації складних функцій за методом Дж. Зойтендейка. При побудові алгоритму застосований 
метод можливих напрямків до вирішення задач чебишевського наближення з додатковими 
обмеженнями. 
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Представлен алгоритм описания овражньтх целевьжх функций, которьюй предлагаєтся использовать для 
детальной визуализации распространения зон поражения сильнодействующих отравляющих веществ 
в случає возникновения техногенной авариий. Такая визуализация необходима и возможна при 
использований геоинформационньх технологий для оперативного реагирования на собьттиє. 
Рассмотренью возможности аппроксимаций сложньх функций по методу Дж. Зойтендейка. При 
созданий алгоритма использован метод возможньх направлений в решений задач чебьшевского 
приближения с дополнительньми ограничениями. 

Ключевье слова: аппроксимация, алгоритм, функция, ограничениє, вектор. 


Не зважаючи на досить великий обсяг робіт з апроксимації складних функцій, 
задача точності таких обрахунків залишається актуальною. Ця актуальність наразі 
проявляється при роботі з програмами, що дозволяють візуалізувати за допомогою 
геоінформаційних технологій зони ураження сильнодіючими отруйними речовинами 
(СДОР). Слід зазначити, що інструменти моделювання дозволяють наносити ці зони 
на карти і схеми у вигляді кола, півкола або сектора, який має кутові розміри 1 
радіус, рівний глибині зараження. Зона фактичного зараження, як правило, має 
форму еліпса включається у зону можливого зараження. Така візуалізація не дає 
картини, що є наближеною до реальності, адже є різні особливості місцевості і хмара 
СДОР не буде чітким еліпсом чи колом. Крім того, у разі хімічних аварій є багато 
речовин, які можуть розповсюджуватися зваженою хмарою, яка заповнюватиме 
певні нерівності території (наприклад, яри) або вкриватиме річкові схили, 
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зависатиме з проміжком від поверхні землі. Такі випадки можна прорахувати, 
враховуючи особливості речовин, напрями розповсюдження |  пересіченість 
місцевості, проте дуже важко представити у вигляді візуальної моделі для осіб, що 
приймають рішення (ОПР). Також візуалізація розповсюдження хмари ЗДОР не 
може бути представлена у вигляді кола чи сектора і у випадку, коли на шляху хмари 
буде річка і хмару СДОР частково потягне за течією. Високі будівлі і споруди на 
шляху СДОР також частково розірвуть контур. 

Огляд існуючих моделей лінійного / програмування, які здебільшого 
використовуються для прорахунків ситуацій, що зазначені, доводить, що ці моделі 
не завжди адекватні реальним ситуаціям. Наприклад, застосування широко 
розповсюджених градієнтних методів може бути неефективним в задачах яружної 
цільової функції, тобто, коли лінії цільової функції сильно витягнуті (мають форму 
еліпсів) у межах оптимальної точки. Такі функції є складними і дослідники 
уникають їх використання -- «дно» таких функцій може бути прямим або звивистим 1 
представляє собою підмножину точок, де поділ на існуючі і неіснуючі змінні зникає, 
тож з будь-якого напрямку функція змінюється повільно. 

Метою роботи є представлення алгоритму опису яружних цільових функцій. 

Слід зазначити, що питання апроксимації функцій свого часу широко були 
досліджені Василенко В.О. (|, Дзядиком В.К. |2|, Поповим Б.О. |3), Люком Ю. |4| 
та іншими українськими і зарубіжними вченими. Проте в роботі використовуються 
елементи підходу, запропоновані Дж. Зойтендейком |5|, що був представлений ще у 
1963 році, проте розглядався дуже мало з причин важкості розрахунків за цим 
методом. У наш час, коли навіть персональні комп'ютери дозволяють вирішувати 
складні задачі, цей метод є цікавим для подальших досліджень та застосування в 
практичній роботі. 

Для представлення алгоритму стисло розглянемо задачу апроксимації деякої 
функції. 

Нехай на проміжку Їа, Б) задана безперервна обмежена функція / (х). 


Розглянемо кусочно-поліноміальну функцію Р(х)є С Ча,5), яка найкращим чином 
5 1 по 
наближує / (х)за підходом Чебишева. Виразом С (а, Б) означимо клас функцій, 


безперервних на відрізку Їа, 5) разом з першою похідною. Явно, що для Р(х) 
матиме місце наступний вираз: 


10) хе а, С | 
509 хе |С,, С, | а) 


Точки а-С, «С, «С, «..«С. «С. -Б приймемо як невідомі. 


Функції РАЕЗ! із-1,5 є поліноміальними зі ступенем не менше 2. Тобто, 
наведена задача, у випадку якщо 70) має однаковий ступінь, і є задачею побудови 


сплайн функції з фіксованими вузлами. 
Задача побудови Р(Х) зводиться до кількох завдань представлення поліномів 


найкращого наближення 70) в розумінні підходу Чебишева до функції / (2) для 


хе Іс, са Щ-ою) Цей факт виходить з принципу оптимальності Беллмана. Саме 


О О.М. Трофимчук, О.О. Кряжич 191 


155 1561 - 5359. Штучний інтелект, 2015, М» 1-2 


тому достатньо розглянути задачу побудови полінома найкращого наближення до 
Рі (х) на деякому інтервалі. Цей поліном повинен ще задовольняти умови, що 


забезпечують відповідну гладкість Р(С). 


Нехай задана функція / (г) і деяка дискретна множина точок: 


вена ТЕ ГоВ| 
У,за, (кції). 


Треба відшукати поліном заданого ступеня К: 
п, суУх ї. 
із0 
який мінімізує величину г(0)стах | у (г, - Пп, (г) по усіх х з області Ле, де 
А визначається: 
да їхев Мапа) 9) пОб) ї-0,1 ) 
іно у«бтчі 


Якщо прийняти, що й за) 


, то задача, що розглядається, буде 


еквівалентною наступній задачі лінійного програмування: 


тіп є 


ах екв) 

щу 

-Уаук-єє- Гб) 
є 

Ха -/(,,) 

Уа о) 

Й 


У "Філ г Гу ) 


ізі 
К 
, ' 
Р 7 діа з 4 (Уа) 
ізі 


є20; віп 


(2) 


Специфіка наведеної задачі лінійного програмування полягає у тому, що 
матриця обмежень А (аз )ї 0,ю; /-0,п--1| має прямокутний вигляд і кількість 


рядків домінує над кількістю стовпців, К «« МУ. Тому для вирішення поставленої 
задачі обрано метод можливих напрямків Дж. Зойтендейка (5|. Через те, що метод 
передбачає наявність нерівності, то умова виразу (2) буде розписана як дві 
нерівності і на майбутнє буде припущено, що всі обмеження (2) мають вигляд 
нерівностей. 

Тепер безпосередньо звернемося до методу можливих напрямів Дж. 
Зойтендейка та представимо алгоритм вирішення поставленої задачі. 

Нехай нам дана довільна задача лінійного програмування: 
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(3) 


Як їі всі методи лінійного програмування, градієнтний метод вимагає 
відшукання точки, яка задовольняє обмеження задачі лінійного програмування. 


Позначимо її Х'- (е ОН ) Тоді для Х 0 виконується: 
Уажльва ху20 і-1, Р: 7 -1,к. (9) 


На відміну від симплексного і двоїстого методів, вирішення задачі лінійного 


програмування Х" може й не бути базисною точкою, що значно спрощує вирішення 
задачі. У цьому дослідженні припустимо, що така точка нам відома. Тоді 


ірраціональна процедура знаходження рішення задачі (3) зводиться до наступного: з 
К 


точки Х" обираємо напрямок 5, за яким величина уа 9; має найбільше значення і 
і 
Кк о 
вектор 5 -(5),..ь5,) задовольняє обмеження У Р.5, «0, ізі Р (РАРУК), де 
уч! 
| | 0 
матриця Ро (Р.) складена з умов матриці обмежень (3), які для точки Х 


виконуються як рівняння, тобто, для матриці Р маємо: 


о за 
2 Вуху зі ізіБ, 
Додаючи сюди і умову невід'ємного невідомого. Після обрання напряму 5, 


обираємо довжину кроку Л для переходу у наступну точку Х/, виходячи з умови, 


що  Х'повинна задовольняти (4). Точка, на якій зупиниться процес, буде 
вирішенням задачі (3). 

Для побудови алгоритму слід більш детально зупинитися на виборі напрямку 
5. Знаходження вектора 5 - (55-55 :.) зводиться до знаходження рішення наступної 


задачі математичного програмування: 


Кк 
4.5. -з тах 
2, Ж (5) 


УР, «0 (-т2) 


до якої, як правило, додають ще одне обмеження (нормалізацію) на вектор 
5 2 (5,58, ). Для дослідження обираємо обмеження: 


К 
нн (7 
узі 


Оскільки застосовуючи запропонований метод на практиці, дуже важко буде 


0 5 5 
вибрати деяку точку Х , яка задовольнятиме (4), то замість задачі (3) можна 
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вирішити задачу, яка у деякому сенсі є еквівалентною задачі (3), тобто застосувати 
метод можливих напрямків до вирішення задач чебишевського наближення з 
додатковими обмеженнями: 


Кк 
У 4 ху - М, тах 


ПІ 


К дог. 

Уазуху Весь ізі, (8) 
іч 

хо 6-20 /2а1К 


де М є великим невід'ємним числом, а величини визначаються системою 


0, якщо Б 20 


і 


-1І, якщо Б. «0, і-1,Р 
Припустімо, значення невідомої 2 дорівнює 2, - тах(- Б/Р «0, із ГР) У 


зазначеному випадку вектор Х б в (00,..0 Є, ) стане початковим вирішенням задачі 


(8). А якщо область умов, що задана у (4), є не порожньою, то задача (3) матиме 


оптимальне рішення, а невідома 2 дорівнюватиме 0. Саме тому, у разі отримання 


від'ємного рішення задачі (8) Х." - їх р оно с зо ми матимемо і оптимальне 


рішення задачі (3) Х,, з ефока очей 5 


Враховуючи, що ентропія інформації описується моделлю, яка визначає 
невизначеність ї повної групи випадкових подій або випадкових станів 
п 
Е 2/-У; р;Їоє р,, а за змістом представляє собою зворотну величину до кількості 
їзі 
інформації, можлива деяка кількість п випадкових подій з імовірністю р)...р,, ЯКІ Не 
відповідатимуть умовам, що прийняті. Наприклад, число М буде числом з 
плаваючою точкою, коли неможливо представити нуль для невідомої 8. 


Враховуючи, що при розрахунках на обчислювальній техніці отримання нуля 
залежить від багатьох факторів, то для побудови початкового вирішення задачі (2) 
можна застосувати додаткові перетворення і розглянути наступну задачу (9). 
Перехід від задачі (2) до задачі (9) обумовлюється тим, що значення точок 7, та 
її (г.), а також обчислення похідних / (а) гі Чь), завжди мають деяку похибку. Саме 


тому, замість рівностей (2), можна обмежитися вимогами виконання відповідних 
умов нерівностей: 


пиа)- (а) « ахє 
' ' 
пиа)- (а) « азє 
Аналогічно будуються умови для точки Р. Величини с,(, додатні та обрані в 
залежності від необхідної точності виконання нерівностей. 
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Уачеє ге) 
Уаок-ава ло) 
ЗУ нано ОО) 
бана -ає Є (у, ) 


і-о (9) 
ї 

5 У амаді п ОСЕСО раРе ) 
і20 


к 
У "Чу лої; -а,є «Г(У,) 
Ра 
к 
7 Ує і -4.ох "ОзЄ2- Го) 
г 


ї 
У "Філмн; зо ГУ) 


ізі 


Е 
- ра і 4уаХ СОС Гу) 
гаю 


є20, із12,.,М 


Вибір початкового вирішення для системи (9) можна отримати наступним 
чином. 
Візьмемо х, «0, і - 0, К . Значення є визначаємо за формулою 


ом ах Оу 0, і-«0О,М-1 
а. 


іє 


де а, - коефіцієнт при є (а,, --1 за, т ). 


У цьому випадку точка Х'-«400,..О є": буде задовольняти обмеженням 
Ен 
задачі (9), тобто, її можна означити як початкову точку. 

Таким чином, вирішуючи задачу (9) методом можливих напрямків, можна 
отримати вирішення задачі найкращого наближення поліномом П С) функції / (2) 
на а, 5) 

Далі припустимо, що поставлену задачу слід вирішити на кінцевому проміжку 
Їж, 81 де функція / (х), що апроксимується, є обмеженою, безперервною та 
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однозначною. При обраному масштабі оберемо відрізок одиничної довжини Ге с) 


Також приймаємо, що всередині цього відрізку відсутні полюсні точки. Попередньо 
досліджуємо на цьому відрізку / (2). Приймаємо, що С,» С,. Нехай відрізок 


розбито на М частин, тобто М -- І точки апроксимації є заданими. Серед них є: 
а) точки С, 1 С); 


6) М фіксованих точок; 
в) (м з 1)- М -2 «М - М - 1 звичайних точок. 
Нехай х; - довільна точка ділення; хо - С, і ху, - С,. Вважаємо, що нам 


. заоі 0 . П 
задані значення / (2) в кожній із точок х,, тобто задана сукупність / (- і У 


50) ГИ ГО), (10) 


о) 
РО НИК МР (1) 


Визначимо 


а також обчислюємо величини 


(0) 


; на відрізку 


Тепер досліджуємо поведінку функції вх), заданої точками с 
ЇС,, С 1 Будемо розрізняти наступні три випадки: 


1) б -с «Є, де 2150 - досить мале число. Тобто, функція сх) З 


відомою точністю поводиться як постійна величина. У цьому випадку природно 
покласти ступінь полінома рівній одиниці п- 1. 


о. Тоді п «2. Якщо ж 


(0) 


є) приймає нульове значення в тт різних точках, між якими перебувають от; 


2) Функція ех) - знакопостійна на с, /ЗА йо 


відмінні від нуля, то п-т 2. 
3) сх) - знакозмінна на Гени 6 Ї З врахуванням того, що число М настільки 
велике щодо довжини відрізка, що ймовірність втрати зміни знака стає дуже малою. 


Нехай вх) змінює знак у / точках з С, о і крім цього в т різних точках 
(0) 


у 70), але не 


сі? - 0 (між цими точками обов'язково перебувають такі, у яких с 


- . - . 0 . 
має знака в їх оточенні. Нехай далі в Р точках, зазначених /чт/ / ( у 0. Тоді 


ступінь полінома 
пзїзтеряі (12) 


Задамося деяким числом є» 0. Будемо говорити, що Р(х) задовольняє по 
точності, якщо 
тах|Р(х)- 705) «Є 
хє а, 5) 
Очевидно, що кожний поліном 70), що становить Р(І), повинен 


задовольняти по точності. 
Нехай на деякому інтервалі с, віє а, 5) побудовано поліном /. (2) ступеня п 


Побудувавши відповідну задачу лінійного програмування, як було показано вище (9) 
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і вирішивши її, ми одержимо величину є. Вирішення задачі відшукання / (х) З 
уточненням ступеня полінома й підінтервалу апроксимації є кінцевою. 

Враховуючи викладене, можна визначити перелік фактичних параметрів, які 
потрібні для використання в алгоритмі, для створення програми опису яружних 
цільових функцій: 

РІ - кількість полюсних точок на РИМА включаючи й кінці інтервалу 

апроксимації. 

Р2 -- кількість фіксованих точок на а, 5. 

МО - кількість точок на відкритому інтервалі (0, Е 1). 

ФІ - кількість пар (х,, /(х, )), які задають таблицю функції /(х) х є |а, 51. 

ЕЇ - довжина підінтервалів апроксимації (попередній крок перебирання 

а, р. 

Якщо виконується протилежна нерівність, то число  підінтервалів 
апроксимації буде більшим, ніж це передбачено вектором полюсних точок С П : РІЇ 
Цим будуть уведені нові полюсні точки. 

02 - задана точність апроксимації (відповідає є ). 

ГР - відповідає величині а, 

Р6 - відповідає величині й, 

С І : РІЇ - масив полюсних точок, записаних у порядку зростання, включає 

точки а Й Б. 
У/ П - РІ) - масив фіксованих точок, записаних у порядку зростання. 
У Й є мо| - масив точок з відкритого інтервалу (0, Е 1). 


ХХ,ЕЕ Й й ФІ| - масиви відповідних точок 

Же а, БТи 7(.,) ла ФІ 

(Сх ре ЕЕ | 1-1... ФІ) 

Власне, алгоритму опису яружних цільових функцій, буде наступним: 

1. Здійснюється введення необхідних параметрів. Лічильники Н, КО, К і 
перший елемент масиву ХХ Іо) представляють первісні значення. 

2. Вибір фактичних параметрів і визначення величин, потрібних для роботи 

програми. 

Визначається величина /) - довжина підінтервалу апроксимації. 

Здійснюється одночасна побудова масивів Х, Е, ЕК, КІ. 


абз 


5. Визначається попередній ступінь 2 полінома апроксимації за 
допомогою масивів Х і К. 

6. Здійснюється зменшення довжини ДО підінтервалу на половину 0: - Р/2 
й відновлення лічильників Н, К, КО. 

7. Побудова матриці обмежень за допомогою масивів Е, 7, КІ. 

8. Побудова СІ й початкове рішення ХХІ задачі лінійного програмування. 

9. Побудова напрямку 51. 

10. Обчислення довжини кроку для побудови нового вектора ХІ. 

11. Побудова нового вектора ХІ. 

12. Блок коректування ступеня. 
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13.Зниження ступеня. Побудова нових обмежень задачі лінійного 
програмування. 

14. Підвищення ступеня. Побудова обмежень задачі лінійного програмування. 

15. Фіксація результатів і перехід до нового підінтервалу апроксимації. 

Таким чином, у наведеній роботі представлений первинний алгоритм опису 
«яружних» цільових функцій, достатній для програмної реалізації мовами 
програмування або за допомогою мови опису алгоритмів 0МІ.. 

Перспективою представленого дослідження є прорахунок конкретних ситуацій, 
які описуються яружними функціями та візуалізація результатів. 

Враховуючи те, що задачі отримання рівномірних наближень сплайнами з 
мінімальною погрішністю розвивалися у багатьох роботах лише у теоретичному 
плані, практична розробка з метою програмної реалізації є актуальною та 
необхідною. Слід зазначити, що наведений підхід та алгоритм може бути 
застосований у сфері підтримки прийняття рішень для вирішення багатьох задач, 
пов'язаних з описом складних об'єктів, розробки програм для пожежних роботів, що 
призначені входити у закриті приміщення для виконання своїх функцій і задач, 
працювати на територіях радіаційного забруднення та виконання інших завдань з 
візуалізації ситуації за допомогою геоінформаційних технологій. 
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ВЕ5СМЕ 
О.М. Тгоїутеспик, О.О. Кгуахрусп 
Дівогіїт дезсгіріїопя сиПу обіесіїує Гипсбіоп5 

Тре аїдогійрта дезсгірйоп5 биПу обіесбїуе Гипсйїопз8. ТРі8 аїєогіфрта сап Бе цп5ед бог 
декаїед уізцайганоп ої Пе дїзітібийоп ої 70пез ої деїеаї ої пієНіу їохіс 5иб5іапсез іп Фе 
еуспі ої іпацзігіа| дізазіег. Тіз гепаегіпє 15 гедиігед мПеп изіпє беоїпіогттайноп 
тесрпоїоєіез Гог гарід гезроп5е їо ап еуепі. 

Тре ц5е ої ргадіепі піеФрод5 пау пої Бе ейесіїує іп Фезсгібіпє Гипсйоп5 єиПу. 
Треве аге Гипсйопя біс Бауе ягеайу еЇопсатед Ппез (еШр5е8) міїбіп Фе Боппдагіез ої 
Фе орійташт роїпі. УУБеп сгеайпо ре аїдогійта ргоро5ед ц5е арргохітайоп ої сотріех 
йипсбіоп5 Бу Ше теїод'я С. 2оціепацК. 
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Тіз рарег ргезепі5 їБе іпіба! аїсогіййт. ТБіє5 аїсогіййт сап Бе іпріетепіса їп 
ргобгаттіпо Іаприаєе5 ог єгарпіса! дезсгірйоп Іапеєцаєє ЮМІ, М/Пеп сгеаїпє Фе 
аїсогіййт иц5ед тефодй роззібіе аїгеспопя їп ргобіет 50Іуїпє СПпербузпеу'я8 арргохітайоп 
гезігіспоп5. 

Саісціайоп ої геа! 5іїшанопя оп їі5 аЇсогіт 15 а ргоплі5іпо, 5(иау. 


О.М. Трофимчук, О.О. Кряжич 
Алгоритм опису яружних цільових функцій 

У статті представлений алгоритм опису яружних цільових функцій, який 
пропонується використовувати для детальної візуалізації розповсюдження зон 
ураження сильнодіючих отруйних речовин у разі виникнення техногенної аварії. 
Така візуалізація є необхідною і можливою при застосуванні геоінформаційних 
технологій для оперативного реагування на подію. 

Застосування широко розповсюджених градієнтних методів може бути 
неефективним в задачах яружної цільової функції, тобто, коли лінії цільової функції 
сильно витягнуті і мають форму еліпсів у межах оптимальної точки. Саме тому, в 
роботі при побудові алгоритму пропонується використання підходу з апроксимації 
складних функцій за методом Дж. Зойтендейка. 

У наведеній роботі представлений первинний алгоритм опису «яружних» 
цільових функцій, достатній для програмної реалізації мовами програмування або за 
допомогою мови опису алгоритмів ОСМІ,. При побудові алгоритму застосований 
метод можливих напрямків до вирішення задач чебишевського наближення з 
додатковими обмеженнями. 

Перспективою представленого дослідження є прорахунок конкретних ситуацій, 
які описуються яружними функціями та візуалізація результатів. 

Надійшла до редакції 01.07.2015 
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